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We compute the magnetic profile of spin Peierls solitons in a simple Heisenberg model with 
magneto elastic couplings, using independently the DMRG method of White and the Hartree Fock 
approximation. Fhe results obtained in such a static model are incompatible with existing NMR 
data on CuGeO, but the distribution of averaged spins agrees qualitatively with the data. We 
conclude that the dynamics in the spin plus lattice system must be included in a more detailed 
theory of the spin Peierls transition in CuGeO. 
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I. EXPERIMENTAL MOTIVATION, MODEL AND TECHNIQUE OF CALCULATION 

Spin Peierls compounds such as CuGeO have solitons in their incommensurate phase |Q| the magnetic profile of 
which can be measured by NMR techniques |2j and which must be accounted for by any theoretical model. Recently 
a systematic study of solitons in the Heisenberg spin Peierls model was done by Quantum Monte Carlo methods [[|, 
with results that differ essentially from those obtained in the xy model Q. However, these authors did not determine 
the distribution of spins and the corresponding NMR signal in the incommensurate phase, which are the subject of 
the present paper. 

Here we calculate the magnetic structure of solitons in the Heisenberg Spin Peierls model via the DMRG method 
of White and also within the Hartree Fock approximation jl2|]. Our (essentially exact) DMRG results in the 
dimerised phase agree with those of ref || and they are also surprisingly close to the results of a Hartree Fock 
calculation. 

The magnetic profiles in the incommensurate phase that we compute by Hartree Fock disagree with those recon- 
structed from NMR data pj. However, if we eliminate the oscillations in these profiles by averaging over even and 



odd sites, we find spin distributions that agree qualitatively with those obtained by NMR in ref[ | 
The spin gap Peierls in CuGeO and similar materials jlj may be due to (i) magnetoelastic coup 
interactions among the chains JO]. We consider the simplest realisation of case (i): 
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lt n are spin 1/2 operators, all energies are measured in units of J and k is the only parameter of this model. The 
AJ„ jn+ i are due to elastic deformations A J n ,n+i ~ (x n — x n +i) the fluctuations of which are thought to be stabilised 
by inter chain interactions and which are treated classically. This approximation ignores the competing energy scales 
of phonons |^] and magnons for CuGeO (an attempt to go beyond the static approximation was undertaken in, for 
example, B). We also ignored second nearest neighbour exchange couplings that might be needed in a realistic model 
of CuGeO l^]. The stationary point of < H > with respect to the deformations AJ n ^ n +i is given by 

AJ n „ tl+ i = < ^„V„ + i > +const (2) 
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and this inhomogenous classical background AJ„ in+1 is the key difficulty in tackling the Spin Peierls Hamiltonian 
of eq(|l|). We use the DMRG method of White in the original spin variables, while for Hartree Fock we use the 
Jordan- Wigner transform |l(| of eq(]l|): 
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As usual, the Hartree Fock (HF) approximation [Q is given by an energy functional 
E\P) = < H spins > 
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plus a consistency condition 
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II. FIRST TESTS OF OUR METHODS 



To test our DMRG |l3| and HF approaches, we consider the ground state energy of eq(Q) for prescribed alternating 
dimerisation AJ at B = T = k = 0. In the HF calculation there is an oscillating effective coupling K n ^ n+ i = 
K + AK(—) n that reflects oscillations in J n ,n+i = 1 + AJ(— )™ and that includes oscillations of the Hartree Fock 
variables. We find (details are given in an appendix) 

K = s + AsAJ, AK = sAJ + As (6) 

As = s = l + -(K-V) 
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where T> ,K, are standard elliptic functions with modulus k = y 1 — {—^} 2 ■ Figure(l) displays, descending in energy, 
(i) the ground state energy of the renormalised xy model 

H xy = + J n+1 ,„(« +1 + slsl +l ) (7) 



(ii) the Hartree- Fock energy and (iii) the DMRG result (DMRG reproduces the Bcthe Ansatz result Eo — 1/4 — log 2 
at AJ = 0). 



ground stale energy as a function of dimerisation 
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FIG. 1. energy per site of a dimerised Heisenberg chain according to, from top to bottom, renormalised xy model, HF 
approximation and DMRG 
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The DMRG and Hartree Fock results are close and at A J = 1 where the system decouples into independent dimers, 
the Hartree Fock approximation is exact. Less satisfactory is the incorrect gap of the Hartree-Fock result at A J = 
that was avoided by Bulaevski |(| by assuming the Hartree Fock parameter to be uniform along the chain. We 
conclude that the Hartree Fock approximation is distinct from the renormalised xy model. 

Both in DMRG and HF we iterate eq(|J) to converge towards the fixed point. We test this procedure on the 
s = 1/2 solitons of the xy spin Peierls model of eq(Q) or, equivalently the charge=l/2 solitons of the Peierls model of 
polyacetylene |{|. We recall that in the xy Spin Peierls model < s z x >= for half the points {{L + l)/2] on a chain of 
odd length L with periodic boundary conditions This is due to the fact that the Hamiltonian couples only even 
with odd points and is of the form 

A 



H xy = 4,+MiP, M = ( " A+ Q ) (8) 



with a spectrum that is symmetric under E — > —E 
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By the symmetry E «-> —E there are 2L eigenvalues of ±E plus the eigenvalue E = of the form 



(9) 
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which alone contributes to < s x > and which is responsible for the vanishing of < s% > on half the lattice. All this is 
borne out by figure 2: 




FIG. 2. spin profile of a soliton in the xy Spin Peierls model 

The magnetic profile of a soliton in the full Heisenberg spin Peierls model as obtained by DMRG and HF are given 
in fig. 3. Clearly, this profile is quite different from that obtained in the xy Spin Peierls model. 
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soliton profile by DMRG and HF 
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FIG. 3. spin profile of a soliton in the Heisenberg Spin Peierls model according to HF and DMRG for ft = 0.75 
The DMRG result gives the symmetric profile in fig. 3 and it is hard to differentiate between the two profiles. 

III. SOLITONS IN THE INCOMMENSURATE PHASE AND THE NMR SIGNAL 



Since DMRG and HF give comparable results for the spin density in the dimerised phase and since DMRG is 
difficult for long inhomogenous chains, we shall use HF to estimate the distribution of spins in the incommensurate 
phase 15 1. 



incommensurate phase, s z and AJ 
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FIG. 4. spins and lattice deformation in the incommensure phase for a chain of L = 200, ^j- = 0.1, k = 1.5 
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We crudely fix the parameter k by requiring that the period of the soliton structure be comparable to that observed 
in the NMR experiments. Consider the configuration of A J and s z , multiplied by an oscillating factor, in the 
incommensurate phase for L = 200, k = 1.5, ^j- = 0.1 in fig. 4. 



distribution of spins in the static model 
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FIG. 5. NMR signal from the static Spin Peierls model 

The NMR data reflects the distribution of spins along the chain or the number of points that carry a prescribed value 
of spin. The configuration of spins as displayed in fig. 4 (modulo an oscillating factor), gives rise to the distribution of 
fig. 5. We see immediately that this distribution is incompatible with the NMR data of ||: it contains positive and 
negative spins and the value of s^ ax — s^ in disagrees by a factor of ~ 10 with the value quoted by 0. 

So we are forced to conclude that the static Heisenberg spin Peierls model for CuGeO is incompatible with the 
NMR data. The most plausible resolution of this 



distribution of averaged spins 
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FIG. 6. NMR signal from the distribution of averaged spins 
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discrepancy is that the motion of a nuclear spin in the NMR experiment averages over the dynamics of the spin 
plus lattice system. To mimic such an average that is, strictly speaking, outside of the domain of validity of our 
calculation, we take averages \ [s z n + s* +1 ) between even and odd spins. The distribution of these average spins that 
is given in fig. 6 is qualitatively similar to the distribution of spins in the xy model, in that the averages have one sign 
only. The value of s^ ax — s^ in for the distribution of averaged spins is of the same order as in the NMR data of [0 , 
but the agreement with the data is only on a crude and qualitative level. 



IV. CONCLUSIONS 



We have calculated the magnetic structure of solitons in a static Heisenberg Spin Peierls model via DMRG and the 
Hartree Fock approximation. Thanks to the fact that the renormalised xy model is not the HF approximation of the 
Heisenbcrgmmodcl we found reasonable agreement between HF and DMRG results in the dimerised phase. 

Within our static model, we obtained a spin distribution in the incommensurate phase that is incompatible with 
the NMR data on CuGeO. We then assumed that the NMR data reflect averages in time over fluctuating spins and 
lattice deformations and mimicked such internal dynamics (that is beyond the scope of our static approach) by taking 
averages over even and odd spins. The resulting spin distribution agrees qualitatively with the NMR data. 

We conclude that the dynamics of the spin plus lattice system must be taken into account in a proper theoretical 
description of the Spin Peierls transition of the CuGeO system. 
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V. APPENDIX 



To correct an error in the litterature, we give the details of the Hartree Fock calculation in the dimerised phase i.e. 
for prescribed and uniform AJ with B = T = =< >= p n ^ n — ~ = 0, see eq(||): 



H eff = E 2 Kn ' n+1 ^+^ n + h -°- 
Kn,n+1 = Jn,n+1 * S n ,n+1, S n ,n+1 = 1 — 2p n ,n+l = 1 — 2 < ?/;„ ?/> n _|-i > 

E(p) = —— Jrun+1 (Sn,n+l s n+l,n ~ 1) 



A = y K 2 cos 2 p + AK 2 sin 2 p 
ip(p) = u(p)a(p) + v*{p)P(p) 
i/j(p + tt) = v(p)a(p) - u*(p)/3(p) 



(11) 
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We are interested in alternating couplings: 

j n , n+1 = i + Aj(-r (12) 

K n ,n+1 = Jn,n+lS n ,n+l = K + AK(— )™ 

AK = sAJ + As, K = s + AJAs 
E(p) = -^volume * (s 2 + As 2 + 2AJsAs - l) 

Fourier transforming the effective Hamiltonian and diagonalising it 

H -ff = \ E ( K + A^R") ri+iipn + h.c. (13) 
= E Hp)[-^ + (p)a(p)+f3+(p)(3(p)] 

\ P \<tt/2 
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we find two consistency conditions: 



P2n,2n+1 = < 1p$ n 1p2n+l >= t + At (14) 
P2n+l,2n+2 = < l/ , 2«+1^2n+2 >= t — At 

t = < T e i "[«*(p)«( P )-t;*(p)»(p)]=-i(K-D) 
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At T = B = both conditions can be combined into one fixed point equation for — 
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1 - 2t, As = -2At (15) 
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From {As, s} we can compute the ground state energy E(AJ). The astonishing quality of the HF calculation raises 
the challenge of going beyond this approximation in a controlled way. 
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